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A theory of thermohydrodynamics in two-dimensional electron systems in quantizing mag- 
netic fields is developed including a nonlinear transport regime. Spatio-temporal variations of 
the electron temperature and the chemical potential in the local equilibrium are described by 
the equations of conservation with the number and thermal-energy flux densities. A model of 
these flux densities due to hopping and drift processes is introduced for a random potential 
varying slowly compared to both the magnetic length and the phase coherence length. The 
flux measured in the standard transport experiment is derived and is used to define a trans- 
port component of the flux density. The equations of conservation can be written in terms of 
the transport component only. As an illustration, the theory is applied to the Ettingshausen 
effect, in which a one-dimensional spatial variation of the electron temperature is produced 
perpendicular to the current. 
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1. Introduction 

Studies on spatial and temporal variations of thermo- 
dynamic and hydrodynamic variables reveal important 
properties of the studied system. 1 ' 2 A length scale and 
a time scale of variations, which occur in response to 
an external perturbation as well as in self-organization, 
are fundamental quantities which are never obtained by 
studying the system in uniform steady states. 

Studies in this direction have already been performed 
by many authors in two-dimensional electron systems 
in quantizing magnetic fields, which exhibit the quan- 
tum Hall effect in the region of low temperatures and 
small currents. 3 ' 4 Theoretical 5 and experimental studies 
on spatial variations of the electrostatic potential in equi- 
librium have revealed a unique spatial pattern consist- 
ing of compressible and incompressible strips, which is a 
manifestation of the Landau quantization. These stud- 
ies have recently been extended to nonequilibrium states 
with the applied current, and the response of the electro- 
static potential to the applied current has been investi- 
gated experimentally 6 and theoretically. 7 However, most 
of these previous works have been restricted to quantities 
related to the chemical potential and the electrostatic po- 
tential. The electron temperature is known to be the key 
variable in the large-current regime where the breakdown 
of the quantum Hall effect 8-10 takes place. The spin den- 
sity is the order parameter of the ferromagnetic phase 
which appears when the Landau level filling factor is an 
odd integer. For a complete description of thermohydro- 
dynamics in quantum Hall systems, studies need to be 
extended to spatial variations of such variables and to 
their temporal variations. 

Spatio-temporal variations of the electron tempera- 
ture have been investigated indirectly by measuring spa- 
tial 11-13 and temporal 9 ' 14, 15 variations of the diago- 
nal component of the resistivity tensor. Kawano and 
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Komiyama 16 ' 17 have recently investigated the electron 
temperature more directly by measuring the local cy- 
lotron emission intensity using a scanning optical probe, 
and have revealed a variety of spatial variations of the 
electron temperature. The first attempt to construct a 
general theory for spatio-temporal variations of the elec- 
tron temperature and the chemical potential in quantum 
Hall systems was made by one of the present authors. 18 
However, this previous theory has been found to employ 
an inexact formula of the thermal flux density. 

In this paper we develop a theory of thermohydrody- 
namics in quantum Hall systems, which is a revision and 
a generalization of the previous theory. 18 We describe 
spatio-temporal variations of the electron temperature 
and the chemical potential in the local equilibrium in- 
cluding the nonlinear transport regime with use of the 
equations of conservation, as in the previous theory. 18 A 
revision is made by writing the equations in terms of the 
total thermal-energy flux density due to the drift motion. 
We generalize the previous theory to study thermohydro- 
dynamics in the presence of a nonuniform potential such 
as a slowly-varying confining potential. 

The organization of the paper is as follows. In §2, we 
introduce a model, assumptions, and macroscopic vari- 
ables. Random potential due to ionized impurities is as- 
sumed to be slowly varying in the scale of both the mag- 
netic length and the phase coherence length. In §3, we 
describe the equations of conservation and introduce a 
model of the number and thermal-energy flux densities 
due to hopping and drift processes. Each of the flux den- 
sities is a linear function of gradients of the electron tem- 
perature, the electrochemical potential, and the electro- 
static potential. The coefficients in these linear functions 
depend on the local electron temperature and the local 
chemical potential. Such dependences are the origin of 
nonlinearity in this theory. In §4, wc derive formulas for 
fluxes measured in the standard transport experiment for 
a narrow sample with a steep confining potential. In §5, 
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we define a transport flux density for each of the mea- 
sured fluxes. We show that the transport flux density 
defined in this manner is equal to the total flux density 
minus the magnetization current density. We find that 
the equations of conservation can be written in terms 
of the transport component only. It is shown that the 
transport flux densities satisfy the Einstein relation and 
the Onsager relation, while the total flux densities do 
not. In §6, as an illustration, we apply our theory to the 
Ettingshausen effect, in which a one-dimensional spatial 
variation of the electron temperature is produced perpen- 
dicular to the current. In §7, we review briefly the exist- 
ing literature on the thermoelectric-transport theory in 
quantizing magnetic fields and on thermohydrodynamics 
in quantum Hall systems. In §8, conclusions are given. 

2. Model, Assumptions and Macroscopic Vari- 
ables 

2. 1 Drift and Hopping Processes 

We consider a two-dimensional electron system (in the 
xy plane) in a perpendicular magnetic field B = (0, 0, B). 
The local potential energy V\ oc consists of the disor- 
der potential due to ionized donors with the fluctuation 
length scale of Z v h ~ 0.1 yum 19 and the macroscopic po- 
tential V . The macroscopic potential V consists of the 
potential in equilibrium such as the confining potential 
and that induced by the applied current. The fluctua- 
tion length scale Z v h is much longer than the magnetic 
length I = \Jhc/e\B\ (e > 0) which is about 0.01 fim at 
|B|=5T. 

An electron loses the phase coherence by a scattering 
from other electrons and phonons. The phase coherence 
length is assumed to be less than l v h- Therefore the An- 
derson localization due to the interference in the disorder 
potential as well as the energy quantization for closed 
orbits around potential hills and valleys arc negligible. 
Then we employ the classical-drift model, 20-22 in which 
each electron state is described by a wave packet with 
extent of the order of I drifting along the equipotential 
line. The energy of the wave packet at r wp in the Landau 
level specified by a set of quantum numbers a — {N, a} 
with the Landau index N = 0, 1, 2, • • • and spin a is given 

by 

el + ViocKp), (1) 

where e° is the energy of the Landau level in the absence 
of the potential. The wave packet follows the equipoten- 
tial line of Vioc with a velocity given by 

dr wp I 2 B / 1\ 

__ = _ Ss eVMoc(r wp ), s B = Wl , e=^_ x Q j. 

(2) 

A wave packet with energy at the center of the broad- 
ened Landau level follows the extended equipotential line 
of the disorder potential and contributes to the macro- 
scopic current by the drift motion given by eq.(2). Other 
wave packets follow closed orbits within each vh-region. 
Here a vh-region is a subspace which is bounded by ex- 
tended equipotential lines and contains one potential hill 
or valley. A hopping of such a localized wave packet be- 
tween neighboring vh-regions by a scattering from other 
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electrons contributes to the macroscopic current. 23 A 
tunneling through a saddle point of the potential also 
produces the macroscopic current, although its contri- 
bution is estimated to be small. 24 Compared to such 
intra-Landau-level transitions, the contribution of intcr- 
Landau-level transitions to the macroscopic current is 
negligible in the slowly-varying potential. 23 In this paper 
we consider the drift and the hopping processes within 
each Landau level in calculating the macroscopic current. 

2.2 Local Equilibrium 

We employ the local-equilibrium approximation, which 
has been used extensively in various systems and is de- 
scribed in detail, for example, in refs.l and 25. We as- 
sume the local equilibrium within each vh-region, which 
means that the energy distribution of electrons in each 
vh-region is described by the Fermi distribution function 
with the electron temperature T c and the electrochemical 
potential /U ec of the vh-region: 

f(e, Mec, T e ) = l/{cxp[(e - n cc )/k B T c ] + 1}. (3) 

Each of T c and /i ec differs from vh-region to vh-region 
in noncquilibrium states. The same electron distribution 
within the condition of constant T c has been used in a 
theory on the quantum Hall effect. 20 To maintain the lo- 
cal equilibrium, it is necessary that the applied current 
and the magnetic field should vary slowly in time com- 
pared with the time necessary for the relaxation to the 
local equilibrium. The slowest process involved in this 
relaxation is the transition of electrons between Landau 
levels in the case where the distribution of electrons to 
different Landau levels deviates from that in the local 
equilibrium. We also assume that phonons are in equi- 
librium with the lattice temperature Tl, and that Tl 
does not change by the presence of the applied current. 

2.3 Spatial Averaging and Macroscopic Variables 
Since the potential due to ionized donors is random, 

quantities such as electron and energy densities fluctu- 
ate in the length scale of Z v h ~ 0.1 fim. We eliminate 
such random fluctuations by a spatial averaging, while 
we retain variations in the scale larger than l v h which are 
produced by the confining potential and by the applied 
current. We assume that T e and fi cc in noncquilibrium 
vary slowly in space compared with Z v h. We note that 
the spatial variation of the measured cyclotron emission 
intensity is also macroscopic because of the limited spa- 
tial resolution of 50 ^m. 17 

Thermohydrodynamics studies spatio-temporal vari- 
ations of macroscopic variables. The macroscopic vari- 
ables in the present system are the electron temperature 
T e (r,t), the electrochemical potential /i cc (r,t), and the 
potential energy V(r,t). The chemical potential fJ,(r,t) 
is defined by fi = fi cc — V. The macroscopic electric field 
E(r,t) is given by E — VV/e. The macroscopic poten- 
tial V(r, t) is determined by electrostatics in terms of the 
averaged electron density which is a function of T c and \x. 
Other macroscopic variables T c and ^ cc are determined 
by the following equations of conservation. 
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3. Formulation with the Total Fluxes 

3. 1 Equations with the Total Flux Densities 

There are two equations of conservation in our model 
of quantum Hall systems. One describes the conservation 
of the electron number, and the other describes the con- 
servation of the energy. The equation of the momentum 
conservation is absent since in the classical-drift model 
the velocity of the wave packet is not an additional de- 
gree of freedom but given in eq.(2) to be a function of its 
position. 

The equation of the electron number conservation is 



dn . 
~0t = Jn ' 



(4) 



where n is the electron density and j n is the number flux 
density. The equation of the energy conservation is 



dt Je 



(5) 



Here e is the energy density and j € is the energy flux 
density, while Pl is the energy loss per unit area at point 
r due to the heat transfer between electrons and phonons 
and is in general a function of jj,, T c , and TL. The energy 
density is the sum of the kinetic energy density and the 
potential energy density: e = ekm + nV. 

Equations describing the time evolution of /x and T e 
are derived from those of n(/j,, T e ) and the entropy density 
s(n, T e ), respectively. The equation for the time evolution 
of s is derived using eqs. (4), (5), and 



T c ds = de — fiocdn, 



(6) 



to be 



ds 



Teg- t =-V- j q - V Mcc -j„-P L , (7) 



where j q is the thermal flux density defined by 

3 q = 3 e ~ MccJn- 



(8) 



3.2 Number and Thermal Flux Densities 
3.2.1 Hopping Components 

First we consider the number flux between neighboring 
vh-regions due to hopping processes, denoted as J^ op . 
Since we consider only hopping processes within each 
Landau level, we have J^ op = J2 a Jn° p - Each Jn° p is in- 
duced by the difference in the electron temperature AT C 
and that in the electrochemical potential A/i cc between 
the two vh-regions. In the first order of AT e and A/i cc , 
4 P is gi ven b y 

Jt p = A a A^ cc + B a AT c . (9) 

In the following we show that the coefficients A a and 
B a are related to each other in our model of hopping 
processes. 

The wave packets contributing to the hopping pro- 
cesses are only those in the vicinity of the boundary 
between the two vh-regions, since the transition rate is 
negligible when the distance between the wave packets is 
much larger than I. The energies of such wave packets are 
confined within an energy range around e a (r) with width 
Thop where e a (r) — £° + V(r) and r hop ~ Fl/l vh with 
T the width of the broadened Landau level. Therefore 



the corresponding occupation probability is well approx- 
imated to be 



fa fi^ai Mcci Te)> 



(10) 



when Thop <C fcB^e- In the case of elastic hopping pro- 
cesses between a vh-region with /j, ec and T c and its neigh- 
bor with ^ cc + A^ oc and T + AT , J^° p is given by 

Jn°a P = -C Q [/(£a,Mcc +A Moc , T c + AT C ) - f(e a , Mcc, T e )]. 

(11) 

In the first order of AT e and A/x ec , we obtain 



(12) 



Also in the case of inelastic electron-electron scatterings, 
this formula is approximately applicable as long as the 
energy change in each hopping process is much smaller 
than k&T e . We assume eq.(12) in the following. 

The hopping number flux density j^ op averaged in the 
macroscopic scale is then given by 



•hop 
J n 



dT c 



(13) 



Here D a is written in terms of the transition rate of each 
hopping process, which depends on the disorder poten- 
tial. The disorder potential in turn is a function of fi and 
T e since it is screened by electrons. 

We introduce transport coefficients V^ x and L^ x and 
write j^ op as 

j h n °? = -Ll\Vn ec - LllT-'VTe. (14) 

From eq.(13) we have 

Lll = e~ 2 a xx = (kaTJ-^Dafail - f a ), (15) 

a 

L x l = (fceTe)- 1 J2 D *fM f*)(e° a m)- (16) 

a 

When we assume the N dependence of D a to be D a = 
(2N+l)D with D a the coefficient for N = 0, the above 
formula of L x \ coincides with that for short-range scat- 
terings in the self-consistent Born approximation in 
the case where T -c k^Tc. Note that the above formula 
is applicable as long as Thop <C k^T c . The coefficients 
L xx and L x xx are functions of T e (r, t) and fi(r,t). In the 
linear-response regime, L xx and L x xx are to be evaluated 
in equilibrium and eq.(14) reduces to the so-called phc- 
nomenological equation. 1 In our model, a nonlinear effect 
is taken into account in dependences of L xx and L x xx on 
deviations of T e (r,t) and t L { r i^) from their equilibrium 
values. 

The thermal flux density [eq.(8)] is given similarly, in 
hopping processes, by 



i hop = 



^2^a - A*oc)Po 



dfg 

d^cc 



dfa 

dT n 



(17) 

This expresses that an electron in the Landau level a 
carries a thermal energy e a — fi cc . In terms of transport 
coefficients, j^ op is expressed as 



•hop 
Jq 



L x lT- x VT c , 



(18) 
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where 



r21 _ rl2 

XX XX ' 



(19) 

Lf x = (fceTe)" 1 £ Q Oa/a(l ~ /«)(<£ m) 2 - (20) 

3. ,2. 1? Drift Components 

The local flux density due to the drift motion fluc- 
tuates spatially since the local potential Vjoc contains 
the random potential. We average the local flux density 
to obtain the macroscopic flux density. The macroscopic 



number flux density in the Landau level a, j 
ten as 



drift 



, is writ- 



adrift 

J na 



(f(s° a + Vioc, Mcc, T c )/i- 1 Ss eVMoc) av • (21) 

Since localized states make no contributions to the 
macroscopic flux density, the occupation probability of 
localized states in the above equation can be replaced by 
that of extended states. Then we have 



•drift 
J na 



!(el + V^ CCl T c )(h- l SB eVV loc ) B 



(22) 



Since the spatial average of VV5 oc is equal to VV, then 



the number flux density j 
becomes 



drift 



due to all the Landau levels 



•drift 



with 



rii _ a yx _ f£ \ " f 

- e 2 - h Z^«- 



Similarly we have for the thermal flux density 

'21; 



•drift 
Jq 



with 



K, 



21 



SB \ -> 
h ^ 



The cocfhcicnts Lli and K?,l can be written as 

y •*> y^ 



2wl 2 



Kf x = s B 



2nf 
h 



-n 



(fi + T e s ) 



(23) 
(24) 

(25) 
(26) 

(27) 
(28) 



in terms of thermodynamic quantities n , fio, an d s o- 
Here fio is the thermodynamic potential density corre- 
sponding to the grand canonical ensemble in the absence 
of disorder (V\ oc = V): 



H (T e ,fj,,B)- 



2ttP ^ 



In 



l+exp(- 



(29) 



while n and so are the electron density and the entropy 
density, respectively, in this case: 

'dQ \ (dn 



n = 



T e ,B 



(30) 



fi,B 



3.2.3 Total Flux Densities 

The total flux densities j n = j h n op + j£ rift and j q = 

jhop + jdrih written as 

j nx = -L x lV x n ec + LllV y V - LllT-'VxTe, (31) 
j ny = -LHV X V - L^V^ec - L^T-^Te, (32) 
j qx = -i^V.Mec + Kf x V v V - LllT-'VxTe, (33) 

jqy = -K 2 y lVxV - L^V^ec - L™T^V v T e . (34) 

The total flux densities are produced not only by V/i C c 
and VT C , but also by VV, and therefore they are in 
general nonzero in equilibrium. The standard transport 
experiment measures a flux through a cross section of the 
sample, which is zero in equilibrium. Such a flux, which 
is induced in nonequilibrium, is calculated in the next 
section. 

4. Measured Fluxes in a Narrow Wire 

4-1 Model and Assumptions 

In this section we derive formulas for the fluxes mea- 
sured in the standard transport experiment. We consider 
a narrow two-dimensional system along the x direction 
and separate the system into the bulk region < y < Ay 
and the edge region y < 0, Ay < y. In the edge region we 
assume the presence of a confining potential which in- 
creases to infinity so that the electron density decreases 
to zero within a length scale of the order of l v ^. The 
width of the edge region An is quite small and varia- 
tions of T c (x,y) and fi ec {x,y) within An are negligible. 
Therefore we calculate the fluxes in the zeroth-order of 
An (zero edge-width model). We also assume that the 
width of the bulk region Ay is small compared with the 
length scale of variations of T c and fi cc , and calculate the 
fluxes in the first-order of Ay (narrow wire model). 

4-2 Edge Currents 

We introduce coordinates (£, 77) for each boundary of 
the two-dimensional system. The unit vector normal to 
the boundary, directed to the outside of the sample, is 
denoted by n. We take the n axis in the direction of n 
and the £ axis along the boundary in the direction of en. 

We calculate fluxes in the edge region n c d E c < n < 
Vcdgc + An. In the region n > r/edge + An the electron 
density and the flux densities are assumed to be negli- 
gible. Since An is small in the present steep confining 
potential, the hopping flux in the edge region is negligi- 
ble. The drift flux is not negligible because of the large 
gradient of the confining potential. With use of eq.(23), 
the drift number flux is given by 

J c „ dgo = K n en, 



(35) 



with 



dn 



11 „, Jv,. 



K n =~-^Y, I dn—f(el + V, fi ec ,T e ). (36) 

Since the r\ dependence of ^ oc and T c is neglected, 
Kn = ^-J2j f(e^cc,T c )de = -s B —n , (37) 
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where e a , /x ec , T c and flo are to be evaluated at rfedge- 
The quantitiy K n is related to the magnetization per unit 
area, M, at r] cdgc by 

#n(<7edge) = (c / e) M (f] cdgc ) , (38) 

since the edge electric current I and M are related 27 by 
J = —cMen. The thermal flux is given by 

J° dgo = K q en, (39) 

with 



S B x f°° 

K q = X] J (e - Moc)/(e, Moc, r e )rfe, 



(40) 



which is written as 

72 



dK n 


- T 11 


dK n 


rl2 

yx 


dfi 


^yxi 


dT e 




dK q 


= K 21 


dK q 


j22 

yx 






dT c 





K q = sb— J (<W) + T cSa (p')) dp'. (41) 

Derivatives of K n and K q are given by 

fits 7\12 

(42) 
(43) 

(44) 
(45) 

(46) 



with 



L 12 



SB h 



T c s = Kyi + K n 



J 22 _ „ 
L vx - S B 



2nP T 2 f 

h e 7_ 



L 22 is also expressed from eqs.(40) and (43) by 



Lf x = 2K q + ^Y.^-^) 2 U- 



4-3 Measured Fluxes 

4-3.1 Fluxes through the Sample 

First we consider fluxes through a cross section of the 
sample along the x direction. Each of the number and 
thermal fluxes is written as the sum of the bulk and edge 
contributions. 

The number flux, denoted as J^ x , is given by 
r&y 

Jl = / 3nxdy + J^ so(A ^ + J^ gc(0) . (47) 
Jo 

Here Jnx 6 °^ and Jf^ &ct " Ay ^ are the edge current flow- 
ing in y < and that in y > Ay, respectively. The 
sum of these edge currents is zero in equilibrium. In 
noncquilibrium states it is induced by the differences 
A/i = (Jb(Ay) - m(0) and AT C = T c (Ay) - T e (0), and 
is written in their first order as 

Jnt ec{Ay) + J^ sc(0) = K n (Ay) K n (0) 

= LiiA M +^AT e . (48) 

The flux density in the bulk region, j nx , is given by 
eq.(31) and is independent of y in the lowest order of 
Ay. Therefore J% x is given by 

Ttr rl2 £12 

(49) 



The measured number flux J„ x is written only in terms 
of V^oc and VT since it is zero in equilibrium. 
The thermal flux J^ x is given by 



rtr ■ j , Todgo(Ai/)tr , 7 cdgc(0)tr 

°qx ~ I Jqxuy T o qx T u qx , 

Jo 



(50) 



where 



edge 
ecJ"nx 



T0dge,tr _ Tcdge _ rodgo _ redgo , / _ \ t, 
°qx — ,J nx H'cc'-'nx ~ ° qx ~ VPoc H'ec)< J i 

(51) 

Here /j,® c is a reference electrochemical potential and lies 
between /i ec (0) and [i cc (Ay). It is the value of the uni- 
form electrochemical potential in a fictitous thin elec- 
trode, which is inserted into the sample to measure the 
fluxes through a cross section at any x. In the first order 
of Ay, the contribution from edges becomes 

r22 

J edge(A,)t r + J odge(0)t r = ^21^ + + 

1 c 

(52) 

where A/z ec = ti ec (Ay) — /U ec (0). The contribution of the 
last term, K n A/j, ec , due to the diamagnetic surface cur- 
rent has been first pointed out by Obraztsov. 28 Then J qx 
is given by 

Ttr T 22 T22 

^ = -L^V^ec +L%V V ^ -^V x T e +^V y T c . 

(53) 

4-3.2 Fluxes through the Boundary 

Next we consider fluxes along the y direction through 
each of boundaries at y = and y — Ay into a fictitious 
electrode outside the sample. 

First we use again the coordinates (£, 77) introduced in 
§4.2 and write fluxes through a segment of the bound- 
ary with length A£ in terms of the bulk and edge cur- 
rents. We integrate eq.(4) over an infinitesimal region 
with width A£ in 7, e dge < V < ffcdgc + A77. The integra- 
tion of the left hand side of eq.(4) is of the order of At? 
and is neglected. By considering additionally the num- 
ber flux through the segment of the boundary, which is 
denoted as AJ*^, we have 



jn V (Ved g e)A^ 



-edge / 



By taking the limit of small A£, we obtain 

Jnr) (?7edge) 



A / tr 



A / tr 



edge 



For the thermal flux, we have from eq.(7) 



AC 



i w (??edgc) - v 4 j; 



edge 



T cdge 



(54) 



(55) 



(56) 



The above formulas are now applied to the boundaries 
of the sample at y = and Ay. At y = Ay, 



-12 



v i J m — v x Iy n — J^yx v xf- ' rp 



V X T 



(57) 



and the measured number flux A J* r „ is obtained to be 



AJ tr 

ny 

Ax 



-Lyl^xHe 



L x L X 2 X 



(58) 
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The same formula is obtained for the boundary at y = 
For the thermal flux at y = Ay, 



V x jf x ^ = Kf x ^ xV L + -^V x T e , (59) 



and the measured thermal flux AJ^ is 



A J, 



<«L- rl2n ,, T 12 X7 it — ^-V T ~ XX Y7 T 

^ x — ^yx VxMec ^ xx Vy/Jec m V X I C - m V y I e , 



r22 
To 



r22 



(60) 



which is the case both at y — and at y — Ay. 

5. Formulation with Transport Fluxes 

5. 1 Transport Flux Densities 

In this section we introduce a transport flux density at 
each point within the large sample by using the formula 
for each of the measured fluxes in the narrow system in 
the previous section. The transport flux densities, de- 
noted as jn(r,t) and j q r (r,t), are defined at each point 
r by 

L 12 



3n(r,t) 



L 1 



(61) 



with 
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(62) 



where L yx = L yx and, since the system is isotropic in 
the xy plane, we have the following symmetry relation: 
L% = L% and L% = -L% with i = 1, 2 and j = 1, 2. 
The right hand side of eq.(61) is equivalent to that of the 
formulas of the measured fluxes in the previous section. 
We note that the coefficients here depend in general on 
x, y, and t through T c (r,t) and [i(r,t). 

Diagonal and off-diagonal transport coefficients have 
different microscopic origins: diagonal transport coeffi- 
cients are due to hopping processes, while off-diagonal 
ones are due to drift processes. Note also that off- 
diagonal transport coefficients are written in terms of 
thermodynamic quantities and do not depend on any 
transition rates. 

It is straightforward to show from eqs. (31), (32), (33), 
(34), and (61) that the transport flux densities and the 
total flux densities are related by 



3 n = 3n + eVM„ , j £ = j* + eVM £ 



(63) 



with 



M n = -K n , M e = -{K q + fx ec K n ) , (64) 

where j* r and j q are related by j q T = jf - /i cc j^, and 
K n and K q are now to be evaluated at a point r within 
the sample. The quantity M n in the above equation is 
proportional to the magnetization M, M n = —(c/e)M, 
from eq.(38). The above equation states that each of the 
total flux densities, j n and j e , is the sum of the trans- 
port flux density and the current density due to spatial 



variations of the magnetization, which is called the mag- 
netization current density. 

5.2 Equations with Transport Flux Densities 

We now rewrite the equations of conservation, eqs. (4) 
and (7), in terms of the transport flux densities. Because 
the divergence of each magnetization current density in 
cq.(63) is zero, the equations can be written in terms of 
Jn and f q T only: 



dn 



v Jn j 



(65) 



ds 



T e - = - V • 3q - V Mcc • ft - P L . (66) 

The term — V/i oc ■ jn is of the second order of the cur- 
rent, while the corresponding term in eq.(7), — V/i c ■ j n , 
contains a first-order contribution due to the drift com- 
ponent in j n when the confining potential is considered. 

The equation for the time evolution of the entropy 
density is obtained from eq.(66) as 

ds 
dt 



v— . i ' ds\ ( ds 



(67) 



where 



7 =T~ 1 7 tr 

J s Jq > 



(68) 



is the entropy flux density, 

f ) = "Jn • T- 1 V^ee - j* ■ T" 2 VT e , (69) 

is the rate of entropy production per unit area at r, and 



is due to the heat transfer to the phonon system. 



(70) 



5.3 Relations between Transport Coefficients 

The Einstein relation tells that the coefficient in front 
of V/i and that in front of in each of the flux den- 
sities are equal. This is satisfied for the transport flux 
densities since they are written in terms of V^ oc = 
V/x + VV. 

The Onsager relation 29,30 tells that 



Lki{B) = Lik(-B), 



(71) 



where L k i (B) is the coefficient in the formula of the kth 
flux density in front of the Ith thermodynamic force. In 
the present case with the entropy production in eq.(69), 
the thermodynamic forces 1 associated with and j q r 
are given by 



2 



(72) 



respectively. Therefore L k i{B) in this case is equal to the 
corresponding coefficient in eq.(62) except the common 
factor T~ x . The coefficients in eq.(62) have the following 
symmetry with respect to the reversal of B. 



L%{T*,»,-B) = L%{T ,p,B), 
L%(T C ,^-B) = -L%(T C ^,B), 



(73) 
(74) 
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with i = 1,2 and j — 1,2. Then we confirm that they 
satisfy the Onsager relation. 

In contrast to the coefficients of the transport flux den- 
sities, those of the total flux densities in eq.(31), (32), 
(33), and (34) do not satisfy the Einstein relation nor the 
Onsager relation. We also note that the coefficient of the 
drift component of the total thermal flux density, K yx , is 
proportional to T c s + f^o instead of T c s [eq.(28)], while 
the corresponding coefficient of the transport flux den- 
sity, Lyl., is proportional to T e s [eq.(44)] and this trans- 
port thermal flux density divided by T e accords with the 
definition of the entropy flux density. 

Smrcka and Stfeda 31 have proved that the transport 
coefficients for number and energy fluxes in quantiz- 
ing magnetic fields are given by integrals of the zero- 
temperature conductivity tensor in the case of elastic 
scatterings (the relation for number and thermal fluxes 
has been written in ref.32). The same relation holds in 
the present case. Each coefficient L 1 ^ with i = 1,2, 
j = 1, 2, a = x, y, and (3 = x,y can be expressed by 

L% = jy^-^ye-^-H^ie), (75) 

in terms of the zero-temperature coefficients 

L x l {a) (e) = I$°\e) = E Q D a 6(e - e a ), (76) 
L#°>( e ) = -LlTie) = ( SB /h)J2 a 8(e-e a ), (77) 
where 8(e) = for e < and 9(e) = 1 for e > 0. 

6. Linear-Response Ettingshausen Effect 

As an illustration, we apply the present theory with 
the transport fluxes to the Ettingshausen effect, in which 
the gradient of the electron temperature is developed in 
the direction perpendicular to the current in the presence 
of the magnetic field. A spatial variation of the electron 
temperature in the Ettingshausen effect in a quantum 
Hall system has been studied previously in the nonlinear 
transport regime. 18 Here we make a linear-response cal- 
culation and investigate spatial variations and quantum 
oscillations of the electron temperature. 

6.1 Model and Equations 

We consider a two-dimensional system in the region 
— W/2 < y < W/2. We employ the zero edge-width 
model in §4.1. The following two cases are studied. In the 
first case the width of the sample W is comparable to or 
larger than the length scale of the electron-temperature 
variation A. The formula of A is given below. In the sec- 
ond case W is much smaller than A, and here only terms 
in the lowest order of W/X are retained. 

We use eqs.(65) and (66) and restrict the calculation 
to the linear-response regime and to steady states. We 
also assume a uniformity along x where flux densities 
and thermodynamic quantities have no dependence on 
x. The exception is jj, ec which has a constant gradient 
along x. The gradient V x [i cc is also independent of y 
since V y S/ x ^ e c = V^V^cc = 0, and is equal to V X V — 
eE x . Then the equations become 

^yJny = , (78) 
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VyjH + Pl=0. (79) 

We employ the simplest model of the energy loss P^: 

PL = C p [T e (y)-T L ] , (80) 

where C p is a constant. The boundary conditions at y = 
±W/2 are 

Jn„ = 0, jfy=0, (81) 

since the fluxes to the outside of the sample are absent 
(AJ* r y = and AJ^ = 0) in cqs.(58) and (60). Using 
eq.(78) and (81) we have 

j*; = 0, -W/2<y<W/2. (82) 

6.2 Spatial Variations 

We substitute V x /i cc — eE x and V x T e = into the 
formula of jg T y , and use eq.(82) to eliminate V y /x ec . Then 
we obtain 

jf y = -A 2 y ieE x - Af y V y T c , (83) 

with 

A yl = L H - LfyiLyly 1 L yl> ( 84 ) 

A vv = \- L m ~ L y\( L y\y lL yV\ T \, ^ ( 85 ) 
The transport coefficients are constant in the bulk region 
in the linear-response regime. By substituting eq.(83) 
into eq.(79), we obtain the equation for T e (y): 

Af y V 2 y T c = C p (T c -n) , (86) 

with the boundary condition: 

jH = -Af x eE x - Af y V y T c = 0,y = ±W/2 . (87) 

This boundary condition immediately shows that the 
electric field E x along the current induces the temper- 
ature gradient V y T . Solving the equation, we find the 
spatial variation of T c to be 

T c (y) - T L = T a [e-(y +w ' 2 V x - e (v-w/2)/x^ (8g) 

The relaxation length of the T e deviation is found to be 

\=(Af y /C v ) 1 ' 2 , (89) 

where A yy > as is proved from the Schwarz inequality. 
The magnitude of the T e deviation is 

T = (\AfJAf y )(l + e- w ' x y l eE x . (90) 

If we assume that D a = (2N + 1)Dq, Tq/E x depends on 
D and also on C p . 

6.3 Quantum Oscillations 

In order to obtain a universal result independent of D$ 
and C p , we consider \7 y T c /V y [i ec in the limit of W <C A. 
In this limit the y dependence of Vj,T and Vj,/i oc are 
negligible from eqs.(88) and (82), respectively. We use 
VxT e = again and eliminate \7 x fj, ec by using j* r = 0. 
We then obtain 

ft, = - B fyVy^ c ~ Bf y V y T c . (91) 

with 

B yy = L yy ~ L yl( L yl)~ lL yl' ( 92 ) 

B fy = \ L fy - L yl( L yl)~ lL yl] T L 1 - ( 93 ) 
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Fig. 1. Electron temperature gradient VyTo perpendicular to the 
current divided by the Hall field V y /j, cc /e (fJ-ec- the electrochem- 
ical potential) is plotted as a function of the chemical potential 
fi for several values of the lattice temperature Tl. 

From the boundary condition j* 1 ,. — 0, we find 

V B T /V v Mec = -Bfyl B ll ■ ( 94 ) 

Here we assume that D a = (2N + l)Do- Then we find 
that V y T /Vj,/i 0C is independent of Dq and C p . The di- 
mensionless quantity fceV ' v T e /Xt ' yfi ec is a universal func- 
tion of /i/hiVc and kBT^/hwc, which is plotted in Fig.l. 
Here the spin splitting is neglected. We find that the sign 
of VyT e exhibits quantum oscillations as a function of /i 
at fcB?L < hw c . Detailed analysis of the Ettingshausen 
effect will be given elsewhere. 

7. Discussion 

7.1 Thermoelectric-Transport Theory 

Nearly half a century ago, the thermoelectric coeffi- 
cients in a magnetic field have been calculated quantum- 
mechanically (the references are cited in some pa- 
pers 28,34,35 ). Several authors encountered a difficulty of 
violating the Einstein relation and the Onsager relation 
when calculating the coefficients for the local current 
density. This problem was due to the confusion of the 
local current density with the averaged current density, 
which is defined as the measured current through a sam- 
ple divided by a cross section of the sample. The Ein- 
stein relation and the Onsager relation should be applied 
to the measured current, or to its averaged density, and 
the local current density is different from the averaged 
current density in a magnetic field due to the presence of 
the diamagnctic surface current (a description is given in 
a textbook 27 and a review paper 35 ). Obraztsov 28,34 has 
shown explicitly that the Onsager relation holds for the 
averaged current density by calculating the total current 
through the sample including the diamagnetic surface 
current in the case of three-dimensional electron systems. 
He has also shown 3 that the thermoelectric coefficients 
in the absence of disorder can be expressed in terms of 
the entropy. 



A current density, which satisfies both the Einstein 
relation and the Onsager relation, can be defined also at 
each point in a sample. It has been defined by j cond = 
j — cV x M and has been called the conduction current 
density, 27,35 the macroscopic current density, 28 or the 
transport current density. 25 Here j is the local current 
density and cVxAf is the magnetization current density 
with the magnetization M . The current density j cond is 
zero in equilibrium by definition. 

Later the discovery of the quantum Hall effect 3,4 has 
stimulated theoretical 25, 33, 36-44 and experimental 45,46 
studies on thermoelectric effects in two-dimensional sys- 
tems in strong magnetic fields. The edge current in a 
two-dimensional system corresponds to the surface cur- 
rent in a three-dimensional system. The thermoelectric 
power has been calculated in a two-dimensional sys- 
tem of noninteracting electrons in a quantizing magnetic 
field in the absence of disorder 36, 37 and in the presence 
of disorder 38, 42 for short-range scatterings in the self- 
consistent Born approximation. 26 Effects of the localiza- 
tion of wave function on the thermoelectric transport 
have been studied 33 within the classical-drift model. 20-22 
Effects of the interaction between electrons on the ther- 
moelectric transport have also been studied. 25 

In this paper we have defined a number transport 
flux density and an energy transport flux density from 
the corresponding measured fluxes in a narrow wire and 
shown that each transport flux density is equal to the 
total flux density minus the magnetization current den- 
sity, namely that corresponding to j — cVxAf . Note that 
the current in this paper is due to only the orbital-center 
motion as in refs.33,34,37, while those in refs.25,27,28,35 
are due to both the orbital-center and relative motions. 
The localization has been taken into account within the 
classical-drift model, 20-22 and the transport coefficients 
obtained in the present paper eq.(62) reduces to those ob- 
tained by Grunwald and Hajdu 33 in the absence of hop- 
ping processes. The transport coefficients in the present 
paper are generalized so that they are applicable to a 
nonlinear transport regime within the local equilibrium 
by including their dependences on the local electron tem- 
perature and the local chemical potential. 

7.2 Hot-Electron Theory and Thermohydrodynamics 

One of the most typical nonlinear effects is the elec- 
tron heating. Uchimura and Uemura 47 have applied the 
hot-electron theory and the self-consistent Born approx- 
imation 26 to explain the electric-field dependence of the 
diagonal conductivity, observed by Kawaji and Wak- 
abayashi, 8 in two-dimensional systems in quantizing 
magnetic fields. 

Shortly after the discovery of the quantum Hall ef- 
fect, the breakdown of the quantum Hall effect has been 
found, 8-10 in which p xx increases by several orders of 
magnitude when the current is increased up to a critical 
value. As a mechanism of the breakdown, a hot-electron 
model, 8,49,50 has been proposed, which combines the 
electron heating and the strong T c dependence of p xx . 
Later the hot-electron model has been supported by the 
observed spatial evolutions of p xx . 11-13 

Gurevich and Mints 49 have proposed a hydrodynamic 
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equation based on the hot-electron model to calculate 
spatio-temporal variations of T c in quantum Hall sys- 
tems. Their equation, however, is restricted to one- 
dimensional variations of T e , and the thermal flux due 
to the drift motion is missing. 

The previous theory by one of the present authors 18 
has proposed a set of hydrodynamic equations for two- 
dimensional variations of T e and [i ec as well as a model of 
flux densities due to drift and hopping processes, which 
is a generalization of the theory for one-dimensional vari- 
ations by the same author. 51 In this theory, 18 however, 
the drift thermal flux density was not correct: the drift 
thermal flux density in eq.(3.13) there was jg Tllt — jq" lt 
where j^ lft is the drift thermal flux density at T c = 0, 
but it should be j^ rift . Equation (3.13) in rcf.18 should 
be replaced by eq.(7) in the present paper and boundary 
conditions cq.(3.14) should be corrected to ■ = 0, 
j* r ■ rib — where rib is the unit vector perpendicular to 
the boundary. 

8. Conclusions 

We have developed a theory of thermohydrodynam- 
ics in quantum Hall systems to study spatio-temporal 
variations in the linear- and nonlinear-transport regime 
in the local-equilibrium approximation. A nonlinear ef- 
fect has been taken into account through dependences of 
the transport coefficients on the local electron temper- 
ature and the local chemical potential. This theory can 
be used to investigate, for example, the electron temper- 
ature distribution in the vicinity of the breakdown of the 
quantum Hall effect in a system with compressible and 
incompressible strips. 

Quantum Hall systems possess several unique features. 
The formation of the Landau levels gives rise to quantum 
oscillations of thermoelectric properties as a function of 
the filling factor. The activation energy for the transport 
plays an essential role in producing the bistability in the 
breakdown of the quantum Hall effect. The large drift 
current is responsible for the observed long relaxation 
length of the electron temperature. 11-13 These features 
will continue to provide a variety of unique thcrmohy- 
drodynamic phenomena in quantum Hall systems. 
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